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We study pure Lovelock vacuum and perfect fluid equations for Kasner-type metrics. These 
equations correspond to a single Ath order Lovelock term in the action in d = 2N + 1, 2N + 2 
dimensions, and they capture the relevant gravitational dynamics when aproaching the big-bang 
singularity within the Lovelock family of theories. Pure Lovelock gravity also bears out the general 
feature that vacuum in the critical odd dimension, d = 2N -|- 1, is kinematic, i.e. we may define 
an analogue Lovelock-Riemann tensor that vanishes in vacuum for d = 2N -\- 1, yet the Riemann 
curvature is non-zero. We completely classify isotropic and vacuum Kasner metrics for this class of 
theories in several isotropy types. The different families can be characterized by means of certain 
higher order 4th rank tensors. We also analyze in detail the space of vacuum solutions for five and 
six dimensional pure Gauss-Bonnet theory. It possesses an interesting and illuminating geometric 
structure and symmetries that carry over to the general case. We also comment on a closely related 
family of exponential solutions and on the possibility of solutions with complex Kasner exponents. 
We show that the latter imply the existence of closed timelike curves in the geometry. 

PACS numbers: 


I. INTRODUCTION 

Lovelock gravity is the most natural extension of general relativity (GR) in dimension higher than four, as it retains 
the basic character of the theory -the equation of motion remains second order despite the action being higher order 
in Riemann curvature. No other purely gravitational theory preserves this crucial feature. Yet another interesting 
feature of GR is the fact that it is kinematic in three dimensions and it turns dynamical in the next even dimension, 
i.e. d = 4. In three dimensions the Riemann curvature tensor can be written in terms of the Ricci so that there exist 
no non-trivial vacuum solution. If we want this property to remain true as we go to higher odd dimensions there is a 
unique choice that corresponds to pure Lovelock gravity mia- The action reduces in this case to a single Nth order 
Lovelock term, with or without cosmological constant [19] in dimensions d = 2N -f 1,2N -f 2. This is the maximal 
order term in the Lovelock series, higher order terms being either topological or zero. There is no sum over lower 
order terms and in particular there is no Einstein term in the action. 

It is possible to define an analogue of the Riemann tensor for iVth order Lovelock gravity, its characterizing 
property being that the trace of its Bianchi derivative vanishes, yielding the corresponding divergence free analogue 
of the Einstein tensor |3]. This is exactly the same as the one obtained from the variation of the iVth order Lovelock 
action term. Eor the appropriate definition of Lovelock-Riemann tensor and zero cosmological constant, any pure 
Lovelock vacuum in odd d = 2N -f 1 dimensions is Lovelock flat, i.e. any vacuum solution of the theory has vanishing 
Lovelock-Riemann tensor mill. Likewise, even for non-zero cosmological constant, the Weyl curvature vanishes in 
3 dimensions, and so does the Lovelock-Weyl tensor for pure Lovelock in all odd d = 2N -\- 1 dimensions. That is, 
pure Lovelock gravity is kinematie relative to the Lovelock analogue of the Riemann tensor in all odd d = 2N -\- 1 
dimensions. 

Kinematicity of pure Lovelock gravity was shown to be true for static spacetimes [2| and conjectured to hold in 
general. This has indeed been proven in general recently mm using a different definition for the Lovelock-Riemann 
tensor due to Kastor |1|- Using purely algebraic properties, it has been shown that this Nth order Lovelock-Riemann 
tensor can be entirely written in terms of the corresponding Lovelock-Ricci in odd d = 2N -|- 1 dimensions. This 
clearly establishes kinematicity of pure Lovelock gravity, as defined above, in all generality. There is no other gravity 
theory satisfying this property. Eurthermore, for kinematicity to hold in all odd dimensions we should restrict the 
pure Lovelock equation to two dimensionalities only, d = 2N -f 1, 2N -\- 2, else this property would be violated. That 
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is, for any given Lovelock order, TV, the pertinent dimensions are only these two. We will be always referring to 
them when mentioning odd and even dimensions in the text. Conversely for a given dimension, d, the order is hxed 
as TV = [(d — l)/2]. Henceforth by pure TVth Lovelock we would refer to this theory in those only two relevant 
dimensionalities, d = 2TV + 1,2TV + 2. 

In this paper we focus on Kasner-type metrics in pure Lovelock gravity. These give a very interesting class of 
cosmological solutions, which are in fact the simplest instances of homogeneous but anisotropic spaces, and have been 
very effectively employed for studying the approach to the big-bang singularity [S]. In the general relativistic case 
in four dimensions, the approach to the singularity oscillates between several Kasner-like phases, contracting along 
two axes and expanding along the third, followed by a switch over to a different expanding direction, and so on. The 
subsequent phases and the transitions between them can be represented as a map from the space of Kasner solutions 
to itself, the Kasner map, that has a very nice geometric structure |20j . In the higher dimensional setting, while 
considering general Lovelock gravity, the leading order behavior close to the singularity is captured by the highest 
order term, thus pure Lovelock gravity would describe the relevant gravitational dynamics in the approach to the 
big-bang singularity within the full Lovelock family of theories. 

The Kasner metrics were also very instrumental in finding the right definition of Lovelock-Riemann tensor verifying 
kinematicity [T] . For pure Lovelock static vacuum spacetimes in d = 2TV -|- 1 both definitions of Lovelock-Riemann 
tensors due to Dadhich [3] and Kastor [4] vanish and the difference became apparent only while studying solutions of 
reduced symmetry. In fact, pure Lovelock Kasner metrics in odd dimensions have zero Lovelock-Riemann in Kastor’s 
formulation whereas Dadhich’s analogue is in general non-zero. 

In the following sections we will analyze pure Lovelock equations for vacuum and perfect fluid spacetimes and obtain 
solutions in the Kasner class. In Sections II and III, we set up the Lovelock framework and write the equation of 
motion for Kasner spaces. Isotropy of spatial stresses is required for both vacuum and perfect fluid solutions. These 
conditions are solved in Section IV, hnding several classes of solutions or isotropy types. Sections V and VI follow with 
a detailed analysis of perfect fluid and vacuum solutions respectively. For each of the vacuum families of solutions, 
we compute higher order curvature tensors that will allow us to actually distinguish between the different families. 
Kastor’s and Dadhich’s Lovelock-Riemann tensors provide an efficient characterization of the solutions. Finally, in 
Section VII we consider a family of exponential solutions, very closely related to the Kasner class, and we conclude 
with a discussion. In the appendix we entertain the possibility of considering Kasner-type solutions with complex 
exponents. Even though in some cases the metric can still be brought to real form, these spaces contain closed timelike 
curves and thus cannot be considered as viable solutions of any gravitational theory. 


II. LOVELOCK LAGRANGIAN AND EQUATION OF MOTION 


The action of Lovelock gravity and the corresponding equation of motion are given by a sum of homogeneous 
polynomial terms in the Riemann curvature, each of them multiplied by a coupling Ck with length dimension 
relative to the Einstein-Hilbert term. Action and equations for these theories are most simply written in terms of 
differential forms 


^ 2^ 

^ ^ {2k)\{d-2k)l 

^ 2k-i 

^^ (2fc)!(d-2fc- 1)!'^ 


Caias-aa A • • • A ^a2k+l gad ^ 


aia2---cid — ic 


i?“ia2 A ... A Jia2k-ia2k ^ gOa^+i /\ . . . /\ gad-l gb ^ 


( 1 ) 

( 2 ) 


In this language, the torsion and curvature forms are defined via Cartan’s structure equations, 

r“ = De^ = de^ + 0;“^ A , (3) 

A , 


for which we have introduced a covariant exterior derivative, D, with the corresponding connection 1-form 0 ;“^, in 
addition to the usual exterior operator d. In order to make contact with the usual tensorial formulation one imposes 
that the torsion is zero and solves for the spin connection in terms of the vielbein. The equation of motion above 
is obtained upon variation with respect to the vielbein only, leaving the spin connection unchanged, as the spin 
connection variation is proportional to the torsion, therefore set to zero. 

Alternatively, one can introduce a set of {2k, 2fc)-rank tensors [4] product of k Riemann tensors, completely anti¬ 
symmetric, both in its upper and lower indices. 


(fe) ]rafcib2'"b2fc _ n[blb2 
^^aia2'"a2fc 


J^2k-lb2k\ 


0102 


02fc-l02fc 


( 4 ) 





3 


With all indices lowered, this tensor is also symmetric under the exchange of both groups of indices, at bi. In a 
similar way we will denote the contractions of K simply as 


{k)mbib2---bj _ 


Uj&ib2-"bjcj+i'"C2fc 
(fe)^Oi02"-OjCj + i'"C2fc 


V J < 2fc . 


( 5 ) 


In terms of these new objects we can now write 

N 

C = — ^ Cfe and 

fc=i 


G\ 





( 6 ) 


In any of the formulations, this reduces to the usual Einstein-Hilbert action (TV = 1) for d = 4 with the first non-trivial 
N = 2 Gauss-Bonnet (GB) correction appearing in five or higher dimensions. Note that the TVth order term in the 
action gives a non-trivial contribution to the equations only for d > 2N -|- 1. The GB contribution to the action can 
be written 


Ggb = - ^Rp.R^'' + R^ (7) 

Varying this term with respect to the metric tensor g^i, (equivalently the vielbein), we obtain 

- 2R%,,RP^ - 2R^^,R% + 

the GB analogue of the Einstein tensor. 

Eor reasons that will become clear later, we will be concerned here with the case of pure Lovelock theories. For any 
given dimensionality d = 2N -|- 1 or d = 2N -|- 2, we will consider the maximal degree iVth order term as the only one 
in the Lovelock series, i.e. cat = I and Cfc = 0 for k ^ N. 

This particular class of theories has some special properties not shared by any other Lovelock theory. Notice that 
it includes Einstein gravity in three and four dimensions. In fact, pure Lovelock theories will preserve some features 
of d = 3,4 general relativity that are not respected by higher dimensional Einstein’s theory [l[6l[2l. As already 
mentioned, in dimension three the Riemann tensor can be written in terms of the Ricci, in such a way that Ricci flat 
solutions are actually completely flat, they have zero Riemann. In dimension four this does not happen and the Weyl 
tensor is in general non-zero. Pure Lovelock gravities as described above generalize this property to any dimension, 
as it is easy to show. In odd d = 2N + 1 dimensions the (2V, 2V)-rank tensor dehned previously can be written 
completely in terms of the Lovelock-Ricci (or the corresponding Einstein), 

mbi---b2N _ ^ ^bi---b2N-\-l , ^a2iv+l /o\ 

^ai...a2N (2N)\ tQj...a2iv + l'-r b2N+l ■ 

which clearly shows that it vanishes when the corresponding Einstein tensor vanishes. However for even dimensionality, 
d = 2N + 2, the above expression is not valid. The analogous would not vanish because it expression involves the (4th 
rank) Lovelock-Riemann, This proves in general, without any reference to any particular solution, that pure 

Lovelock gravity is kinematic in all odd d = 2N 1 dimensions. In case of non-zero cosmological constant, it can be 
easily verified that the Lovelock-Weyl tensor, the traceless part of vanishes as is the case for the usual Weyl in 

three dimensional Einstein gravity. 

In all relevant odd or even dimensions the whole information about the {2N, 2A^)-rank tensor is contained in 
so we will be referring to this quantity only for all practical purposes. We will be comparing the results obtained for 
this Lovelock-Riemann tensor with those of other 4th rank tensors defined below. 

In particular for TV = 2 we may define two 4th rank tensors (see [1] for details and definitions for any N) 

= R^''°‘^ Rp^^-y + + RR^^-y (9) 

and 

^ ^ (JO) 

Related to the first of these objects, one of us [3] defined an alternative Lovelock-Riemann analogue 


{2)-paP —{2)pap _ 

%J "VV - V'^ 


( 11 ) 
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where TZ is the complete contraction of and c? = 5, 6. Kastor’s Lovelock-Riemann [3], in turn, can be simply 

written as 

^ ^ Q (12) 

In Sec. VI, it would be shown that for a particular class of pure GB vacuum solutions in odd d = 5, would 

be be non-zero while would be actually vanishing. 


III. KASNER METRICS 

Consider a d-dimensional metric with flat t = const, slices, for which each of the n = d — 1 spatial directions scales 
differently with time in a polynomial fashion. 


— dt^ 


(13) 


This is known as the Kasner metric and is a homogeneous but anisotropic space. These type of metrics were also 
considered long ago by Deruelle [8] in the GB case (see also 0EQ] for a more general account of anisotropic models 
in Lovelock theories). 

For this class of metrics, each additional power of the Riemann curvature adds a factor t~^ to the Lovelock- 
Riemann tensor and the strength of singularity increases accordingly with the Lovelock order. Focusing in the 
dynamics near the big-bang singularity located at t = 0, this means that the dominant contribution will come from 
the highest order Lovelock term in the action. Therefore, pure Lovelock theories capture the leading order behavior 
in this regime. In fact, Kasner metrics, as written above, will not, in general, be exact solutions of the equations 
of motion when considering terms of different curvature order in the action. If we want to analyze the subleading 
behaviour of the metric we may keep the next lower order Lovelock term and expand the metric further with two 
polynomial terms in time, gu ^ -I- . Lower Lovelock terms become relevant for sufficiently long times, i.e. 

dd,N-k/C-N ~ CN-kt^^ ~ 1 - 

Since for pure Lovelock the relevant dimensions are only two d = 2N + 1, 2N -|- 2 we would therefore be considering 
the corresponding cases, n = 2N and n = 2N + 1, respectively for odd and even dimensions. In what follo ws, we write 
the Lovelock-Einstein tensor, corresponding to the Vth order Lovelock term, for the metric (13) for odd and 

even dimensions. Just the diagonal components are non-zero. 

For odd d = 2N + 1 dimensionality, there are 2N exponents pi and the spatial components of the Lovelock-Einstein 
tensor are 


_ N{2N-2)l 

^ * f2N 


2N - 1 - Pm 

\ m—1 





(14) 


whereas the time component yields 


Q{N)t ^ _ 


N{2N-1)\ 


Pi-" P2N ■ 


(15) 


Note that does not involve pp, i.e there are only 2N — 1 exponents Pj^a, in the last product or the sum. 

In particular, for d = 5 GB these expressions reduce to 


and 


\ ^ (3 - P2 - P3 - P4 )P2P3P4 


= -^PlP2P3P4- 


In even d = 2N + 2 dimensions the expressions are similar, yet a bit more involved, 

^ / n. \ ( 2., , . 

PyPi] +P]---pi\ j <k-- - <l (16) 




2N-1 


2N 








5 


The number of pj in each product term in the sum is 2N — 1, and there are ( 2 ^^i) = 2A^ such terms. Besides, the 
time component reads 


Q(N)t ^ _ 


N{2N -l)\ 



■Pi = 


N{2N - 1)! 

w 


2 N +1 


-pi ■ ■ ■P2N+1 


y 1 


j < k ■ ■ ■ < I 


(17) 


2 N 


for which each term in the sum is a product of 2N exponents pi. Each product involves a combination of all the 
2N + 1 Pi, i.e there are {^ 2 N^) ~ ^ terms. 

For GB in d = 6 we may write 


^ [(3 - P 2 - P 3 - P 4 - Pb){P 2 P 3 .Pi + P 2 P 3 P 5 + P 2 P 4 P 5 + PsPiPb) + P 2 P 3 .PiPb] 

and 

W2)t 12. X 12 /I 1 1 1 1 

^ t = - -^(PiP 2 P 3 P 4 + P 1 P 2 P 3 P 5 +P 1 P 2 P 4 P 5 +P 1 P 3 P 4 P 5 ■^P 2 P 3 P 4 P 5 ) = - T^PlP 2 P 3 P 4 P 5 -^-^-^-f- 

t t \Pi P 2 P 3 Pa P 5 


IV. ISOTROPY CONDITIONS 

Even though Kasner spacetimes are anisotropic in general, in many situations we are interested in the isotropic 
case for which all the spatial components of the Lovelock-Einstein tensor are equal. This is the case for vacuum or 
perfect fluid solutions that will be the focus of our analysis in the following sections. We shall therefore analyze the 
isotropy conditions 

GWv _ = 0 (18) 

for every pair i,j. Not all equations written in this way are independent, only n — 1 are, and we can choose, for 
instance, a set with fixed i or j, say i = 1 . 


A. Odd dimension d = 2N + 1 


For odd dimensionality these conditions read 


' 2 N 


Q(N)r^ _ QiN)j_ ^ I _ 2iV + 1 j Pi ■ -J Pm^ =0, 

\fc=l / 2 N -2 

In particular fr five dimensional pure GB, it would reduce to 

G(2)1i _ G'(2)22 ^ _P2) (p^ +P2 +P3 +P4 - 3) P3P4 = 0 

and two similar equations for and — G^'^^'^ 4 . 


(19) 


( 20 ) 


B. Even dimension d = 2N + 2 


In turn, in even dimensions we find 


/ 2 Af+l \ 

G^n)k - ~ (p, - p^.) f ^ p, - 2V + 1 j ^ 


Pi"-^ = 0, 

2 N -2 


( 21 ) 


where each product in the sum contains 2N — 2 different pi that are a combination of all {2N — 1) available exponents 
with I 7 ^ i,j, i.e. there are = 27V — 1 of those terms, 

in six dimensions, N = 2, we shall write 

(^(2)11 _ ( 7 ( 2)22 ^ (p^ _P2) (P4 +P2 +P3 +P4 +P5 - 3) (^ 3^4 + papg +P 4 P 5 ) = 0 (22) 

and three similar equations for — G ^^^^4 and . 
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C. Isotropy types 

There are several ways of solving the isotropy conditions. In both odd and even dimensions these reduce to a 
product of three factors. Therefore, for each pair i,j at least one of those factors must vanish. Depending on the 
choice of the vanishing factors we will have several families of solutions or isotropy types. 

• The first family -type (a)- corresponds to the trivial solution pi = pj for every i,j, i.e. all the exponents are 
equal and the metric trivially isotropic. 

• Another very simple solution -type (b)- corresponds to putting to zero the second factor, that in fact does not 
depend on the pair i,j that we pick. 


d-l 

2N-l-Y,Pk=0 ■ (23) 

fc=i 

A single condition ensures isotropy in this case. Whereas (a) type solutions form a one parameter family of 
metrics, (b) type spaces form a n — 1 parameter family of solutions. 

• If we choose to set the third factor of the isotropy conditions to zero instead -type (c)-, for every pair i 7^ j we 
get different conditions in odd and even dimensions. In odd d = 2N + 1 we have to solve 

Pi---Pm = 0 I - ■ -m ^ i,j , (24) 

2N-2 


whereas in even d = 2N + 2 dimensions the condition is 



2N-2 2N-1 


(25) 


The form of the solutions is similar though. As we will see below, in odd d = 2N + 1 this istropy type implies 
that two exponents vanish, say pi = P 2 = 0, whereas in even d = 2N + 2 three exponents have to vanish, 
Pi = P 2 = P 3 = 0 for instance. In both cases, these solutions correspond to vacuum spacetimes and will be 
considered in detail in Section VI. In the N = 1 case with d = 3,4, 2N — 2 = 0, and hence these conditions 
become vacuous. There are no solutions in this class for Einstein-Hilbert gravity. 


• The last possibility -type (d)- is to combine isotropy types (a) and (c) in such a way that for some pairs i,j 
we have pi = pj^ whereas for others the third factor in the corresponding isotropy condition is the one that 
vanishes. As we will see below, this will only be possible in even d = 2N + 2 dimensions. 


Isotropy types (a) and (b) are simple enough and do not require any more explanation. Let us discuss type (c) 
and (d) in more detail. In odd dimensions enforcing just one of type (c) conditions readily implies that one of the 
exponents vanishes, say pi = 0. Automatically all the components of the Lovelock-Einstein tensor vanish except for 
G^i, but isotropy implies that the remaining component must also vanish, and hence another Pi/i = 0. The rest of 


the exponents may take any value, 
which (241 reds, say for j = 1, 


To be more concrete let us focus for a moment in d = 5 pure GB gravity, for 


P2P3 = 0 , P2Pa = 0 , PsPi = 0 


This implies that two of the pi exponents with i 1 must be zero. All the components of the Lovelock-Einstein tensor 
vanish in that case, therefore all solutions in this isotropy type are vacuum spacetimes. This is also true in general, 
that is, isotropy type (c) in all odd d = 2N +1 dimensions implies vacuum spacetime. We will comment more on this 
later on. 

The even dimensional counterpart of the above statement is a bit more complicated to get. For d = 6 pure GB 
again, we have for i = 1, 


PiP5 + P3P5 + P 3 PA = 0, PAPb + P2PA + P2Pb = 0, P2P3 + P2P5 + P3P5 = 0, P3PA + P2PA + P2P3 = 0 


This implies that three of pi must be zero. Note that, if we set one of pi to zero, say, pi = 0, then the isotropy 
conditions would read as 


P2iPAP5 + P3P5 + P 3 Pa) = 0, PsiPAPS + P2PA + P2P5) = 0, Pa{P 2 P 3 + P2P5 + P3P5) = 0, Pb[P 3 PA + P2PA + P2P3) = 0 . 
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up to a factor ^Pi — 2N + 1 that we asume to be non-zero. Otherwise the solution would already be considered in 
the type (b) class. From the previous equations it is now obvious that two more p^’s are zero. Again, this would lead 
to vacuum solutions which will be analyzed further in the corresponding section. 

For general N, imposing the above condition (25) just for one pair i,j implies either two pi^ij exponents are zero, 
for instance pi = P 2 = 0 , or else 


1 

— = 0 

Pk 

In the first case, pi = P 2 = 0, all components of but yield zero. Furthermore, isotropy requires the 

rest must also vanish which implies one more Pi^i ^2 = 0 and we have again a vacuum solution. The other possibility 
X) ^/Pk^i,j = 0, tacitly assumes that none of the exponents involved is zero. However, when considering another 

pair such that pi ^ Pj,Pi, the only way G\ — = 0 is that at least one Pk^i,j,i = 0, thus contradicting the initial 

assumption. 

The only remaining possibility is that all p^’s are organized in two groups such that either p^ = pi 2 with pi ^ p 2 . 
This corresponds to isotropy type (d) and we just have to impose one restriction 

- G^^ ~ ^ 1 = 0 . (26) 

fe^l .2 

The rest of the isotropy conditions are either equivalent to this one or trivial. For pi ^2 having multiplicities ni .2 with 
ni+n 2 = 2N + 1, this reduces to 



rii — 1 no — I 

-H- 

Pi P2 


= 0 


(27) 


with ni ^2 7 ^ 1- For d = 6 GB gravity, we have for instance ni = 3, 712 = 2 and pi = — 2 p 2 , and it is easy to check that 
the isotropy condition G^^ = G\ is verified. 


V. PERFECT FLUID 

Once the isotropy conditions are imposed the equation of motion coupled to matter has the form = 
diag(G*(, GL,..., GL), which matches precisely the form of the stress-energy tensor of a perfect fluid with energy 
density p = —G*( and pressure P = G\. Moreover, all components of the Lovelock-Einstein tensor G\ scale as 
with time, and consequently the same behavior applies to the fluid p and P supporting these geometries. Energy 
and pressure are, in this way, proportional to each other and satisfy a linear barotropic equation of state of the form, 
P = uj p with w = w(pi; d, N). 


(a) All equal pi exponents 


We have density and pressure as given by 


= -G^S = 


(d- 1)! /pWAf 


())■ 


and 




where d = 2N + 1, 2N + 2 respectively for odd and even dimensions. The equation of state is 

f 2N I \ 


which gives for p = 2 {d^i) ^ 5 ^ dust, radiation and stiff fluid respectively. It is an FRW flat model. 










(b) EtlP* = 2iV-l 


Note that the . component is free ofpi. We can however substitute in Eqs (14) and (16) 2A^—1 —= Pi 

so that in both cases has the same form of except for a numerical factor. Thus, upon substitution, all 

spatial stresses become equal, and it turns out that 


P = 


1 


2N -V 


for both odd and even dimensions. 


(c) d — 2N + 1 zero exponents. 

As explained before this isotropy type corresponds to two and three zero exponents in odd and even dimensions 
respectively. This, as already mentioned, are vacuum solutions and thus have vanishing energy density and pressure, 
p = P = 0. 


(d) Pi = pi _2 with multiplicities ni 2 

For the isotropy condition (c), this is the only case of non-vacuum solution and this class of perfect fluid solutions 
exist only for even d = 2N -|- 2 > 6 dimensions. As before we have to further impose the condition, 


Pi P2 

with ni ,2 1 and ni + n 2 = + 1. Substituting this into the equations of motion we get 

G\ -(2iV- +P 2 ) (28) 

whereas the spatial components are 

G\ = G\ ^ pT~"pT~" {2N-1- (m - l)pi - (n2 - 1 )^ 2 ) (29) 

Thus we get an equation of state, 

/ 2A^ - 1 - (ni - l)pi - (712 - 1)P2 \ ^ 1 - Ay 

V {2N - l){pi + P 2 ) E Ay ^ ■ 

The nice parametrization above corresponds to pi = (ni — l)x, P2 = — (^-2 — l)x and A = rii — 71 - 2 . Notice that this 
class of solutions allows very large values for pi ^2 (equivalently x) in which case we obtain P = —p. Conversely for 
small pi ^2 — > 0 we have an almost vacuum solution with P « ). The parameter w = P/p thus spans the whole 

range [—l,oo) in this case. In principle this parametrization would allow for fluids with lower barotropic parameter 
w but these would violate all of the usual energy conditions. Notice that the barotropic constant depends on two 
parameters in this case, one discrete, ni, and the other continuous, pi. For each possible value of 1 < ni < 2N (or 
A) we have a continuous function, w = u}{pi). 


VI. VACUUM SOLUTIONS 

Vacuum spacetimes can be considered as a subset of the more general perfect fluid class of solutions from previous 
section. Sincein all cases we have a equation of the type P = ujp, we just need to impose one more condition, p = 0, 
to get the vacuum solutions. Clearly a perfect fluid or vacuum solution for pure Lovelock would not in general satisfy 
the isotropy condition for the Finstein tensor. 

In odd d = 2N + 1 dimensions p = 0 implies that at least one of the exponents has to be zero. This fact combined 
with the isotropy conditions yields two types of solutions. Considering the isotropy type (b) we get solutions of the 
form 

2N 

Pi = 0 and P 2 = ‘2N — 1 — ^ pi . 

i=3 


( 30 ) 










9 


We have already mentioned that all isotropy type (c) solutions are vacuum. These have not one but two zero 
exponents, say 

Pi =P2 = 0. (31) 

In accordance with the discussion above, these solutions will be referred to as type (b) and (c) vacuum metrics 
respectively. In the way defined so far, these two families have a non-empty intersection. To be precise we will 
consider as type (b) metrics those with only one vanishing exponent, pi>i ^ 0 , in order to have mutually exclusive 
categories. 

These are all non-trivial vacuum Kasner solutions of the theory in odd dimension. Type (a) solutions are trivial 
as all the exponents are equal and thus p = 0 implies that the metric is just Minkowski. Moreover, there are no type 
(d) vacuum metrics in odd dimensions. In fact, there exists no solution with all non-vanishing exponents pi in odd 
dimensions. Besides, in the type (b) case at least one of the exponents has to be positive whereas the sign of the 
rest is a priori unconstrained. There are no further restrictions on the exponents, Pi>2, of type (c) solutions, on their 
signs or otherwise. Both vacuum families of solutions have the same number of free parameters, we have 2N — 2 free 
exponents, Pi>2- 

For even d = 2N + 2 dimensions, there are two ways of solving the vacuum condition p = 0. We may have two flat 
directions, say pi = P 2 = 0 , or 

2N+1 

— = 0 . 

Pt 

This can be combined with the different isotropy types to give three different families of solutions. The first two families 
are analogous to those found in the odd dimensional case and correspond to having at least two zero exponents. For 
type (b) metrics we have [ 21 ] 

2N 

Pi=P 2 = ^ and P 3 = 2A^ - 1 - (32) 



whereas type (c) solutions are all vacuum as explained before. These have three flat directions instead of just two, 

Pi=P 2 =P 3 =0 (33) 

These two families of solutions will be referred to as vacuum types (b.l) and (c) respectively. We again further 
constraint type (b.l) metrics to have only two flat directions, i.e. Pi >2 7 ^ 0, so that there is no intersection between 
types (b.l) and (c). 

For the last family of solutions, corresponding to the case of all pi being non-zero -only compatible with type (b)-, 
we have 


2Af-|-l 

Pi = 2N — 1 and 

i=l 


2Ar-i-i 


-r± 

E - = o- 


(34) 


This family will be referred as (b.2) and is the familiar Kasner metric in its most generic form. Clearly, there have to 
be both positive and negative exponents pi, at least one of each kind. For type (b.l), as it happens in odd dimensions, 
at least one exponent has to be positive but there may be none of negative sign. There are no constraints on the signs 
for type (c) metrics. 

Again type (a) solutions are trivial and correspond to Minkowski flat space and, even though we have type (d) 
metrics in even dimensions, p — 0 would require either pi 2 = 0 or pi -|-p 2 = 0, both conditions incompatible with Eq. 
(27). The latter condition would imply ni = n 2 , but this is impossible as the sum of the multiplicities has to be odd, 
ni + n 2 = 2N -|- 1. Therefore, there are no vacuum type (d) solutions. All this discussion is summarized in Table |l] 

In odd d = 2N + 1 dimensions we have two {2N — 2)-parameter families of solutions whereas for even d = 2N + 2 
dimensions we have three sectors, two {2N — 2) and one {2N — l)-parameter families of solutions. It is clear that 
for d = 3,4 type (c) vacuum solutions are just the trivial flat solution while for d > 4, we have non-trivial vacuum 
solutions in that class. 


A. Curvature tensors 

We may now compute the different tensorial quantities defined in the Introduction, namely Kastor’s Kabcd and the 
alternative TZabcd-, but also the Riemann tensor. Rated', for the different vacuum types. For simplicity we will denote 



10 


type 

d 

isotropy cond. 

vacuum 

vac. curvature 

(a) 

any d 

Pi ^ p , Vi = 1, 2 ... d — 1 

p — 0 

A = 7^ = K = 0 

(b) 

-b.l- 

-b. 2 - 

2N + 1 

2N + 2 

2N + 2 


Pi = 0 

Pi = P2 = 0 

y^d-l J_ _ Q 
Z^i=l Pi 

R, 7?., R 7 ^ 0 

(c) 

2N + 1 

2N + 2 

Pi = P2 = 0 

Pi = P2 = P3 = 0 

all 

7^ = R = 0 

(d) 

2N + 2 

Pi ~ Pi ,2 with multiplicities ni ,2 
"1-1 -h ""-1 = 0 , ni+n2=2N+l 

PI P2 , L : 2 

none 



TABLE I: Classification of isotropic and vacuum solutions. (#i) 7?. = 0 as well for the subset of type (b) ((b.l) in d = 2A + 2) 
with all non-zero exponents equal to one, Pijti ,2 ~ 1- (# 2 ) For flat Kasner (naively type (c)) we also have R = Q. for N = 2, 
except for these exceptional cases, all the tensors that are not in the table are non-vanishing. 


these 4th rank tensors just as K, TZ and R respectively, not to be confused with respective contractions. Note that 
Lovelock vacuum is defined by Rat = R-ab = 0. 

The computation of these tensors is very simple for type (c) vacuum metrics for which both tensors are actually 
zero. This is easy to understand as in both tensors there is at least one set of antisymmetrized 2N indices. For type 
(c) metrics, for both odd and even dimensions, we have at most 2N — 2 non-zero exponents, i.e. 2N — 1 non-flat 
directions including time. As a consequence, at least one of the 2N antisymmetrized indices will have to be on one 
of the flat directions. As any component of the Riemann tensor involving that direction is zero, for all dimensions 
and Lovelock orders we have TZ — 0 and = 0 (also then Ai — 0). Another way of stating the same thing is that 
the spacetime has reduced effective dimension 2N — 1 (number of non-flat directions), with d — 2N -|- 1 flat directions 
added to it, and the corresponding Lovelock-Riemann is therefore zero. The Lovelock-Riemann tensor is non-trivial 
only in deff > 2A^. It would have to be d^s > 2A^ -|- 1 for the Lovelock-Einstein to be non-zero. Moreover, as diseased 
in [T] , in deff = 2the Lovelock-Riemann is completely determined by the Lovelock scalar, whereas for deff = 2N + 1 
it can be given in terms of the Lovelock-Einstein. Note that the Lovelock scalar is proportional to the trace of the 
Lovelock-Einstein except when deff = 2N. In that case the Lovelock-Einstein vanishes but the corresponding scalar 
in general does not. It does vanish in vacuum as we will see below. Vacuum spacetimes have non-trivial Lovelock- 
Riemann only in dimension, des > 2N + I, below that threshold all vacua are Lovelock-flat; i.e. the corresponding 
Lovelock Riemann is zero. 

Eor type (b) solutions in odd dimensions (or (b.l) in even dimensions), a similar argument holds. The effective 
dimension is a priori des = 2N in this case. In principle we have enough indices so that the antisymmetrization 
does not necessarily yield zero as before. We may realize however that, if all exponents are either zero or one, all 
components of the Riemann tensor with any time index vanish, = Pi{pi — l)t~^ = 0, thus effectively removing 

that direction as well. We reduce again deff = 2N — 1, thus implying 7^ = K = 0. Notice that all such solutions belong 
to types (b.l) and (c) (depending on the number of zeros and ones) in both the relevant odd and even dimensions. 
One can easily check that, at least for d = 5,6 pure GB, type (c) metrics, and those of type (b.l) ((b) in 5d) with 
all non-zero exponents equal to one, are the only solutions verifying TZ — 0. 

Something similar happens for for the Lovelock-Riemann. In this case, as the effective dimension is deff = 2N, this 
tensor can be written completely in terms of a scalar, basically the Lovelock term for that effective dimension [T]. 
Also the equation of motion will be given in terms of this invariant, therefore it has to vanish. We can explicitly check 
that by computing the effective Lovelock scalar. 



(27V)! 

{2N - l)t^^ 


' 2 N -1 

Y, P^-^N + l 

_ 


2iV-l 


n Pi 


(35) 


where we have set to zero, Pi> 2 Ar-i = 0. For type (b) metrics this quantity is zero, thus, the whole Lovelock-Riemann 
tensor is also zero. Note that the Lovelock-Enstein is not trivial in this case as the real dimension is bigger than 2N. All 
of type (c) and type (b.l) metrics are Lovelock flat solutions. This had to be the case in odd d = 2N -\-1 dimensions 
because of kinematicity, but it is a non-trivial statement in even dimensions. Contrary to TZ the Lovelock-Riemann 
tensor vanishes for all type (b.l) solutions not just for those with = 1- 

Summarizing, the least restrictive condition among the ones we use is Lovelock flatness. The Lovelock-Riemann 
tensor, M = 0, singles out type (b.l) (all type (b) in odd dimensions) and (c) metrics. 7^ = 0 is verified just for all 
of type (c) but just a subsector of type (b.l) with all non-zero exponents equal to one. These spaces effectively have 
reduced effective dimensionality, deff < 2N — 1, thus any tensor with 2N antisymmetrized indices vanish for these 
metrics. 
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p^ 



FIG. 1: Lovelock flat solutions in GB gravity in the P 4 = 0 plane in d = 5 (also in d = 6 with ps = 0). The vertical and horizontal 
planes correspond to each of the exponents being zero, Pi = 0 for i = 1,2,3, whereas the remaining one is pi + P 2 + Ps = 3. 


Moreover, the Riemann tensor vanishes only when all exponents are zero except for one that may be zero or one. 
Both possibilities correspond to flat spacetime, the former being just Minkowski space in its canonical form. The latter 
-so called flat Kasner—, even though naively belonging to a different isotropy type, also corresponds to Minkowski 
space, a patch of it, in a different set of coordinates. Therefore the previous solutions are, in general Lovelock flat but 
not Riemann flat. Besides, only type (b.2) vacuum solutions have a non-trivial Lovelock-Riemann tensor. One can in 
fact verify that, the Lovelock flatness condition for 6-dimensional pure GB metrics precisely implies that these have 
to have at least two flat directions, thus belong to types (b.l) or (c). The other tensors are also non-vanishing for 
(b.2) metrics. We can use the previous results concerning the vanishing of the different curvature tensors to actually 
classify all vacuum Kasner solutions into the corresponding families. For that we just have to be careful in identifying 
the exceptional cases, i.e. type (b) metrics with all non-zero exponents equal to one and flat Kasner. 


B. Structure of vacuum solution space: the Kasner Shamrock 

The space of vacuum solutions can be easily visualized in the space parametrized by the exponents. Except for 
type (b.2) that has a more complicated form, all of type (b) and (c) vacuum solutions correspond to intersections 
of planes in (we will denote this space K^”^), two in odd d = 2N + 1 dimensions, and three in even d = 2N + 2 
dimensions. These are the Lovelock flat solutions. For the simplest case, N = 2, GB gravity in d = 5,6 dimensions 
we can visualize this in the three dimensional parameter space spanned by {pi,P 2 ,P 3 ) with the remaining exponents 
set to zero. We have in this way four families of solutions corresponding to pi = 0, p 2 = 0, pa = 0 -type (c)- and 
Pi + P 2 + P 3 = 3 -type (b) or (b.l) in odd and even dimensions respectively- (see Figurej^. 

The remaining case, type (b.2), has a much richer structure and will be explored in detail in what follows. These 
solutions verify the isotropy condition (b), 2N — 1 — P* = which implies that the volume of t = const, slices 

grows as This is also true for vacuum type (b.l) spaces. In addition, type (b.2) solutions are the only ones 

that have a non-trivial Lovelock Riemann tensor, they are not Lovelock flat and are in this sense more generic. In 
fact they are in the most generic form of Kasner, with all non-vanishing exponents. 

Type (b.2) vacuum solutions verify one less condition than the other two families present in even dimensions. In 
fact, we will see that type (b.l) solutions appear at the boundary, in some appropriate sense, of the bigger (b.2) 
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family. When taking two of the exponents to zero, the second of (341 conditions have to be understood as a limit 


and, depending on how this is taken, we may solve the condition for any value of the remaining pi. This condition is 
thus vacuous and we recover the (b.l) type solutions. These will appear as codimension one limiting surfaces in the 
space of (b.2) solutions. Notice that type (c) metrics are completely different in this sense, they cannot be recovered 
in this way as they correspond to different isotropy types. All type (b) metrics lie in a plane 'Y^Pi = 2N — 1 in 
whereas type (c) metrics do not. 

As we have two constraints on the exponents, we can always solve for two of them in terms of the remaining ones. 
Defining two new parameters 


2 N 


2 N +1 


1 


c = Pi + P 2 = 27V - 1 - Vpfe and ^ = — + — = - V — 


k =3 


Pi P2 


we can determine pi and p 2 as solutions of a quadratic equation 

p^-Cp+7=o 


yielding 


Pi = 


1±a/1- 




P2 = 


iTi/l- 


ce 


all the remaining exponents Pk^i ^2 being free. There is only one limitation to the above formula, the product of the 
two parameters cannot be in the (0,4) range, therefore 


Ce>4 or CC<0, 


(36) 


otherwise the values the pi _2 become complex. This constrains possible values the remaining pi can take. For a pair 
of complex conjugate exponents we can make a complex change of variables making the metric real. These are also 
solutions to the pure Lovelock (also Einstein) equation, however they contain closed timelike curves and thus do not 
give physically interesting well defined spacetimes. We will comment more on this in the appendix [j 

Another interesting limiting case of the above formulae is ^ = 0 which necessarily implies ^ = 0 as well, unless 
Pi = P 2 = 0 in which case we recover a type (b.l) family of solutions. For C = ^ = 0, we can still choose the limiting 
value of C/'x/C in such a way that we obtain, pi = —p 2 = p, a one-parameter family of solutions. 

Alternatively, we may try to see what the space of allowed P 3 ,P 4 , • • ■ looks like for fixed ^ and C, similar to how the 
Kasner solutions of Einstein gravity can be visualized as the intersection of a plane and a sphere 

^p, = 1 ; ^P? = l- 

In three dimensions these two constraints reduce the space of nontrivial solutions to just one, one of the pi parameters 
being zero and the other equal to one. This is what is known as the flat Kasner solution, as it has zero Riemann 
curvature. In four dimensions we can split the solutions into those with two vanishing exponents, the other being one 
-type (b.l), again flat Kasner in this case- and those with all pi nonzero -type (b.2)- for which the sphere constraint 
is equivalent to the above Y ^/Pi = Oj with i = 1, 2, 3, when the equation for the plane is used. The only type (c) 
metric is trivial in this case. The reality conditions (36) in this case simply amount to the exponents, all three, being 
in the range Pi G [—1/3,1]. The upper bound is obvious from the sphere condition, otherwise at least one of the other 
exponents has to be complex. One obvious difference between d = 4 Einstein gravity and higher dimensional pure 
Lovelock is that the exponents are bounded as ||pi|| < 1 in the former case, equality meaning that the solution is flat 
Kasner. This seems to have some relation to the stability of these solutions m- 

In d = 6 GB gravity we can visualize the three dimensional space of P 3 ,P 4 ,P 5 in a similar way with equations. 


P3 + P4 + P5 = 3 - C 


—I-1— — —C- 

P3 P3 P3 


For even dimensions higher than 6 -and corresponding Lovelock order N > we can always decompose the space of 
exponents pi into lower dimensional subspaces similar to the one above. This will allow to suitably represent the space 
of solutions and also to exploit the rich symmetry structure of these equations. We now analyze the six dimensional 
case in detail. 

To make the analogy closer to the d = 4 Einstein case we can rescale the exponents, qi = Pij (3 — Cfl), such that 


<73 + 94 + 95 = 1 


- + - + - = e(c-3) 

93 94 95 
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0 1 2 3 4 5 


FIG. 2: Level diagram for the function K{X,Y) = S'® defined in the text. The orbits of solutions for the 

different values of k correspond to the level sets of the function K, i.e. K{X,Y) = k. The numbers on the color scale 
correspond to ^ arctan(fc). 


The case C = 3 has to be treated separately. Notice that for ^ = 0 we recover precisely the d = 4 space of solutions, 
the Kasner sphere -a circle in this case-, whereas in general it gets deformed. To see how the deformation parameter, 

= C(C - 3) = (— + —) (pi +P2 - 3) , 

\Pi P2 ) 

affects the space of solution we can plot these spaces. In order to simplify the visualization we project onto the plane 
d 3 + 94 + <75 = 1 parametrizing 53 = ^ = i+x+V^Y ^ _ i- 2 X terms of the new {X,Y) 

coordinates (see Figure [^. Cold colors correspond to negative values of k whereas warm colors indicate positive k. 
We see the thick line corresponding to the Kasner circle {k = 0) and the white lines to one of the exponents being zero, 
gi = 0, or A: — ± 00 . These lines correspond precisely to type (b.l) metrics. For k to be infinite we need ^ —>■ ±00 
(we cannot take C —> cxd otherwise the rescaling would not be well defined) and thus either pi or p 2 vanish as well. 
The intersections of the white lines correspond to two qi = 0, the remaining (b.l) solutions, for any value of k in this 
case (equivalently any value of ^ 1 , 2 )- We can, in this way, analyze the whole set of type (b) -both (b.l) and (b.2)- 
vacuum even dimensional solutions. To complete the analysis one just has to include type (c) solutions that cannot 
be represented in this way. 

Another interesting set of solutions is A: = 1, for which the orbit again yields three lines in the space of solutions. 
These correspond to each of the exponents being one, = 1 for i = 3,4,5, that automatically solve the constraints. 
This same triangle plays a prominent role in d = 4 Einstein gravity as the vertices can be used to define the Kasner 
map. This map is an application of the Kasner circle on itself and corresponds to the evolution of type II Bianchi 
models that have Kasner asymptotics both to the future and the past. The iteration of this map represents the so 
called Mixmaxter attractor, the subsequent transitions between Kasner epochs as we approach the big-bang singularity 
in more general (Bianchi IX) models (see [H] for a recent discussion). It is reasonable to think that a similar map 
may exist for Lovelock gravity models as well, despite the fact that the previous chaotic behaviour disappears in 
higher dimensions for vacuum Einstein gravity [HHIs]. The characterization of the space of solutions performed in 
this section can thus be taken as the starting point for the more general analysis of big-bang singularities for this class 
of theories. 

Interestingly enough, the previous representation makes apparent the existence of three axes of symmetry that 
correspond to the exchange of a pair of exponents, qi -O- g_,-. We have three possible such exchanges. We have also 
another symmetry that corresponds to a rotation of 27r/3 that in turn represents a cyclic permutation of the three 
exponents. Using these symmetries we can always restrict the values of X and Y to one of the six fundamental 
domains, each corresponding to a different permutation of (pa, _P 4 , _P 5 ). 

The representation of Figure is useful, for instance, to discuss the signs of the exponents. In the center of the 
inner triangle all three qi are positive, whereas each time one crosses one of the white lines one of the exponents 
changes sign. We can see thus that over the Kasner circle we always have one negative and two positive exponents. 
This is also true for any solution corresponding to a positive value of k, except for the piece of the orbit inside the 
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inner triangle (only for k > 9) that has all three exponents bigger than zero. Negative k orbits split in two parts as 
well. The piece contained inside the Kasner circle k = 0 still has two positive and one negative exponents, whereas 
the pieces in the outer blue triangles have two negative directions. 

The only solutions that cannot be represented in the {X,Y) plane as above correspond precisely to ^ = 3. In this 
case the system reduces to 


P3+P4+P5 = 0 ; —^-^— = -C 

P3 Pi P5 

and can be treated in a very similar manner. We can again rescale the exponents, qi = —^Pi, so that now we set to 
one the parameter in the second equation, 

93+94 + 95 = 0 ; V—-—h — = 1 (37) 

93 94 95 

Using a similar projection as above, qs = 94 = q^ = now on the plane X]fc 3 9* = 0, we can 

represent this orbit as shown in Figure]^ below (in blue). Besides, if we change the signs of all the exponents qt we get 
the corresponding curves in the complementary regions (in red). Notice that the same system can be obtained from the 
generic case rescaling the g’s (or equivalently X, Y) by k and taking the limit fc —> 0. These curves represent somehow 
the boundary of the {X,Y) space above, brought to finite distance. Remark that in the original representation the 
9i = Pi/{3 ~ C) diverge as ^ 3. Besides, we can also view this as a Kasner sphere with inverted coordinates, 

qi 1/ qi- 

Instead of just for ^ = 3, we could have used the representation in terms of qt and the corresponding rescaling from 
the beginning. We would then had gotten, 

111, 

93 + 94 + 95 — ^ j V———h — — 1 

93 94 95 

where the value of fc = ^(C — 3) is the same as before. In principle, any point can be represented in two equivalent 
ways related as qi = kqi. In addition to this, we can also get one system of equations from the other through an 
inversion, -f-)- I /qi. This inversion can be then rephrased as a symmetry of the original equations, whether in terms 



FIG. 3: Representation of the Kasner circle upon inversion as defined in Eq. (371 with 53 


^ _ X-V3Y - _ 


Qi 


95 = 


-2X 

3 • 
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of Qi, Qi or Pi- This symmetry amounts to 


P- o — 

c - 3 ip. 


or 


Qi ^ -r~ 
kqi 


or 


Qi 


k 

Qz 


(38) 


We can find similar symmetries for any subset of the exponents except for the complete set. 

For any non-zero value of fc, this symmetry would generate from one piece of the orbit some other piece with the 
same number of positive and negative exponents (for positive k) or all opposite signs (for negative k). The only points 
for which we cannot choose the representation we use are those with k = 0. Depending on whether ^ = 3 or ^ = 0, 
the corresponding orbit can just be plotted using qi or qi respectively. Still we can easily obtain one of these orbits 
from the other through the inversion, qi o l/qt- 

The representation in terms of the qi is otherwise not very useful. The first representation, as already seen, is 
projected on the plane 93 + 94 + 95 = 1 whereas for the second we can project on 93 + 94 + 95 = k. Notice, that 
contrary to the original visualization scheme, now each point in the (X, Y) plane does not have a unique value of k 
associated to it, but three. There are three orbits with different values of k through every point. It is much more 
complicated to plot and not very enlightening. 

Remarkably enough, the inversion symmetry qi ^ l/{kqi), when translated into the {X,Y) plane, admits a very 
clear and elegant geometric realization. To see this, instead of performing just the inversion, we will define three 
different inversion maps, each of them being the composition of the inversion with one of the three possible exchange 
symmetries, Qi -f-)- qj with i ^ j. Remember that these correspond to reflection symmetry about each of the three 
reflection axes of Figure We will refer to them as fi with i = 1,2, 3; e.g. fi corresponds to the composition of the 
inversion with the 92 ^ 93 exchange. The other maps are defined analogously. 

The composition of these reflection symmetries with the inversion, each of these maps, verify a very interesting 
geometric property. The straight line on the {X, Y) plane connecting the original point with its image via the map 
always passes through one of the intersections of two white lines, 9 ^ = 0 or fc —±00 (see Figure]^ in green). 
Conversely, the intersection of one such straight line with the orbit of given k value would give me two points related 
by the map. This is very reminiscent of the way one can define geometrically the Kasner map, the only difference is 
the points used to trace the lines. The Kasner map has the intersections of the fc = 1 lines as focal points, instead 
of those oi k —>■ ± 00 . If I pick a different green point to trace the line the resulting points would be related through 
a different one of the maps. Moreover, from a given point, tracing the lines through the three green points I would 
get the corresponding three images. These images are all image under inversion of the same point, up to a reflection 
symmetry. Therefore they are related through the composition of two such reflections, i.e. they are related by a 27r/3 
rotation. Graphically they form an equilateral triangle around the origin of the {X, Y) plane. 

Using this map we can restrict to the values of {X,Y) that are inequivalent under the inversion. This will allow 
us to focus on a compact set of points in that plane, despite the original space of solutions being non-compact. The 
boundary of that compact set (shaded region in Figure]^ is given by the fixed points of any of the maps. It is easy 
to verify that one such set of fixed points corresponds to the circle, 

{X -if +Y^ = 1 

whereas the rest can be obtained by rotating this set by ±27r/3 (see Figure]^ bottom right). Each of the regions 
contained in the shaded region generates under inversion one of the corresponding regions outside. More precisely, 
given the set of solutions for a given k inside the three circles we get its counterpart outside. We can actually perform 
this operation in a completely geometric way just tracing lines. Given three points related by 27r/3 rotations inside the 
shaded region, one can trace the nine lines through these and the green points mentioned before. This nine lines will 
cross in groups of three outside the shaded region generating the images of the three original points under inversion 
(see Figure]^. These are also related by 27r/3 rotations. In this way, given the constant k orbits in the shaded region 
we get the corresponding orbits outside. 

This operation is well defined for all points on the {X,Y) plane, except for the points on the Kasner circle that 
get mapped to infinity (the lines are parallel in groups of three), corresponding to the orbit C = 3, and points on 
the white lines. Each of this lines gets mapped to a single points, one of the green points used to generate the map. 
Notice that these are (b.l) solutions that do not enjoy the inversion symmetry. 

The more involved part of the shaded region corresponds to the inside of the inner triangle. In that region the circles 
overlap dividing the triangle in six regions, three narrow and three wide regions. The map thus maps the narrow 
regions into the wide ones and viceversa. We have chosen the wide regions to belong to our shaded fundamental 
domain but we could as well have chosen the other ones. The shape of this fundamental domain consists in three 
equal leaves of almost circular shape, joined at the center, reminiscent of a shamrock. We may call it the Kasner 
shamrock. Each of the leaves of the shamrock can be divided in two halves, each of these halves contained in one of 
the six sectors from which one can reconstruct the space of solutions with reflections. Thus, from just one of these 
half-leaves we can reconstruct the whole space of solutions using reflections and inversions. 
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FIG. 4: Graphical representation of the inversion maps fi. Any point and its image correspond to the intersection points of 
straight lines through each of the green points and a given fe-orbit. The blue lines correspond to the orbits of solutions with 
the same value of k. For /s the related green point corresponds to qi = q 2 ~ 0. An example of such map appears in the larger 
figure. The rest of the figures show all three images under the inversion maps for points in different regions of the solution space. 
The boundary of the Kasner shamrock corresponds to the points that are mapped to themselves under any of the inversion 
maps (see bottom right corner figure). 


VII. EXPONENTIAL SOLUTIONS 

Another class of solutions closely related to Kasner spacetimes is that of exponential type metrics of the form 

n 

ds^ = —dt^ + ^ e^^'*dx^ 

i=l 


( 39 ) 
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FIG. 5: Graphical representation of the inversion maps of three points related by ±27r/3 rotations in the {X,Y) plane (on 
the left). The tracing of the nine straight lines joining these points (in red) with the points generating the maps (in green) 
determine the three corresponding image points (in orange). The only points that do not have an image in this way are those 
in the Kasner circle whose lines are parallel in groups of three (on the right). 


These exponential solutions are a generalization of de Sitter spaces, with different Huhhle parameters in every direction. 
This is very similar to the Kasner metric in a different form except that the dt^ lacks the customary in front. In 
fact, both Kasner and exponential solutions can be treated together using a slightly more general form of the metric, 

n 

ds^ = -e^^otdt^ + ^ e^^'*dxj , (40) 

i=l 

that reduces to the Kasner metric for Hq = 1 and to the exponential form for Hq = 0. For more general values of 
Hq we can always make a change of coordinates of the form t = (accompanied by a rescaling of the spatial 

coordinates) to bring the metric again to the Kasner form with exponents pi = Hi/H q. The exponential solutions 
can thus be pictured as living in the boundary of the space of Kasner solutions. We can take a limit where pi —> oo, 
at least for some of the Kasner exponents with, Pi/pj = Hi/Hj. 

The role of the Hubble exponents. Hi, is very similar to the pi of previous sections. There are just two differences. 
First while the pi are dimensionless quantities the Hi have dimension of inverse length. This is the reason the (2, 2) 
Riemann components scale as for Kasner whereas they are constant in the exponential case. Second, the way 
these constants enter the Riemann curvature is just slightly different, the t Hf while the remaining 

components, R^^ij = —> HiHj, are equal in form without the scaling factor. For the more general ansatz 

the form of the ti components would be Hi{Hi — Hq) making the connection between the different parametrizations 
obvious. 

At the level of the pure Lovelock equations of motion, we realize that these become of homogeneous degree in the 
new exponents, Hi. From this basic observation and the discussion of Kasner solutions we can readily get the general 
form of the constraints for the exponential metrics. The algebraic form of the equations for the Hi parameters is the 
same as for the pi when we keep just the highest degree in all expressions. Lower order in p terms correspond to 
terms with Hq factors in the exponential representation that are then set to zero. For instance we can take the above 
classification of vacuum solutions for pure Lovelock gravity and write: 

In odd d = 2N + 1 dimensions we have two {2N — 2)-parameter families of solutions 

u Hi = H 2 = 0 -type (c). 

• Hi = 0 and H 2 = — J 2 'i ^3 Hi -type (b), with i7i>i ^ 0, so that there is no intersection between types (c) and 

(b). 

whereas in even d = 2N + 2 dimensions we have three sectors, two {2N — 2) and one {2N — l)-parameter families of 
solutions. 
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• Hi = H 2 = Hs = 0 -type (c). 

» Hi = H 2 = 0 and H 3 = - Yl'i= 4 ^ -type (b.l), again with i?i >2 ^ 0 . 

• Hi ^0, \/i with H^=0 and 3^ = ^ “type (b.2). 

This is summarized in Table We can do the same for perfect fluid solutions or any other. Notice that this class 
of vacuum solutions does not exist in Einstein gravity without cosmological constant. Revisiting the Kasner circle 
adapted to this case we get 


Y.H, = 0 ; = 0 ( 41 ) 

i i 

and the second condition readily implies that all the exponents are zero. 

In order to go to the boundary of the space of solutions, as explained above, we may take pi _2 00 , therefore, 
either ^ 0 or () —)• 00 . In the latter case we can rescale qt = —pilC (Pi —> 00 ) to get the equation of the k —)■ ±00 

orbit 


d-l 

9 * = 1 

The second possibility is to rescale then qi = —pi^ 


d-l 


and ^ 

Qi 

2=3 


(pi —>■ 00 ) to get 


(42) 


d-l 


d-l 


5Z 9i = ?C ^ 0 and 


i=3 


2=3 


(43) 


This is the image of the Kasner circle under inversion qi —)■ l/qi (see Figure]^ for the d = 6 case). Moreover, taking 
both limits at a time, ^ 0 and C —?> 00 we can actually fix the value —>■ /c getting the corresponding orbit, 
analogous to the Kasner case, and its image under inversion. The structure of the space of solutions is very similar to 
the Kasner case. Notice that the inversion Hi 1/Hi is also a symmetry of the exponential space of solutions for the 
whole set of exponents Hi in this case. For other subsets we can define an analogous transformation similar to that 
of previous sections. We can actually identify the parameters qi with a rescaled version of the Hubble parameters Hi. 

We can readily analyze the form of the curvature 2-form (equivalently the Riemann tensor) for the more general 
ansatz (40) above. 


= H,{H, - Ho)e 


A e* 


= Hi Hi 




(44) 


and notice that for every value Hq 7 ^ 0 the Riemann components diverge either to the past or to the future (at the 
singularity) in which case the dominant term in the Lovelock equation will correspond again to pure Lovelock. The 
exponential case, Hq = 0, is quite different as the components of the Riemann tensor become constant and all Lovelock 
terms contribute to the same order in the equation of motion. Besides, every curvature scalar will be constant and 
there is no curvature singularity at all. The geometry is completely regular. 

In the case of Einstein gravity we can add a positive cosmological constant to get solutions of this type. In the 
general case we would get a complicated set of polynomial constraints involving all of the Lovelock couplings. Even 
in the pure Lovelock case, the addition of just a cosmological constant term would complicate the equations to a high 
degree. 

In analogy with the Kasner metrics we can also discuss the values the different curvature tensors take for vacuum 
exponential solutions. We summarize the situation in Table [h] The only difference is that the subset of sector (b) for 
which 7Z — 0 disappears. This is easy to understand as it would correspond to Hi = Hq in the general parametrization, 
thus Hi = 0 and they reduce to type (c) metrics. For this class of solutions we have a perfect hierarchy of solutions 
with regard to which curvature tensors vanish. This is true for N = 2 but expected to hold in general. Type (a) 
vacua have R — TZ = R = 0, type (c) 7^ = K = 0, i? yf 0, type (b.l) (or (b) for d=2N-|-l) has M = 0, i?,7^ yf 0 and 
finally for type (b.2) all those tensors are non-zero. The Riemann tensor is only zero for the trivial Minkowski metric 
in this case -there is no analogue of flat Kasner. We can then distinguish the different families of vacuum solutions 
just by using these tensors! 

In the case of Kasner metrics this is not true due to the existence of exceptions; i.e. flat Kasner, that belongs to 
type (c) yet has R = Q, and the subset of type (b) metrics with all non-zero exponents equal to one, that have 7^ = 0. 
Once this cases are taken care of separately, the above classification scheme carries over to Kasner as well. 
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type 

d 

isotropy cond. 

vacuum 

vac. curvature 

(a) 

any d 

Hi = H, Vi = l,2...d-1 

o 

II 

in 

R = 72 = K = 0 

(b) 

-b.l- 

-b.2- 

2N + 1 

2N + 2 

2N + 2 

Et/ = 0 

Hi =0 

Hi = Hi = 0 

Z^i^l Hi ^ ^ 

R = 0, it, 72 7^ 0 

R = 0, i?, 72 A 0 
i?, 72, R 7^ 0 

(c) 

2N + 1 

2N + 2 

Hi = Hi = 0 

Hi=H2=H3=0 

all 

7^=:M = 0, i?7^0 

(d) 

2N + 2 

Hi = ill ,2 with multiplicities ni ,2 

^ + UfU = 0 , ni + n2 = 2A + 1 

none 



TABLE II: Classification of isotropic and vacuum exponential type solutions in pure Lovelock gravity. 


VIII. DISCUSSION 

Kasner solutions play a fundamental role in the analysis of big bang singularities. In this note we have analyzed 
and classified Kasner type metrics in pure Lovelock theories of gravity. Regardless of the phenomenological interest 
of such theories, they capture the leading order dynamics in the approach to the singularity within the Lovelock class 
of theories. Close to the big-bang the curvature is generically diverging, thus the leading contribution will be that of 
the highest order curvature term, thus pure Nth order Lovelock gravity in d = 2N + 1, 2N + 2. 

Analyzing the conditions for isotropy we were able to classify perfect fluid and vacuum solutions in several families 
denoted as (a)-(d) (just (a)-(c) in the vacuum case, with subtypes (b.l) and (b.2) in even dimensions). In vacuum 
the different types differ mainly in the number of flat directions, pi = 0, the metric has. Type (a) vacuum metrics 
have all directions flat, it is just Minkowski. Type (c) metrics have at least d — 2N + 1 flat directions, thus two and 
three in odd and even dimensions respectively, whereas type (b) solutions have at most d — 2N vanishing pi. In odd 
dimensions these have one flat direction but in even dimensions they may have either two or none, corresponding 
to subtypes (b.l) and (b.2). Lovelock non-flat vacua emerge only in even d = 2N + 2 dimensions only when all 
exponents piS are non-zero, these satisfying the conditions '^Pi = 2N — 1 and '^p^^ =0. This is true for any pure 
Lovelock theory in vacuum. 

Parallely, we analyzed the values the different Lovelock-Riemann analogue tensors take for the different vacuum 
types. We found a nice correspondence between the different families of solutions and the vanishing of particular sets 
of these 4th rank tensors, in a hierarchical way. Type (a) solutions verify R — TZ = M. = 0, type (c) have 72. = K = 0 
and R 0 (except for flat Kasner that has 7? = 0), while type (b.l) -or just (b) in odd dimensions- yields K = 0 
and i?,72 ^ 0 (except for the solutions with all non-zero exponents equal to unity). Finally type (b.2) solutions are 
the ones for which all these tensors are non-vanishing. We can then use this tensors to classify our solutions in the 
different families, taking into account the exceptions. This may be very helpful in classifying solutions in case these 
are not given in the canonical Kasner form, but in some other possibly complicated set of coordinates. Besides, this 
classification scheme may be relevant for more general classes of solutions. 

A similar classification scheme exists for another class of solutions closely related to Kasner. Exponential type 
metrics as those studied in section VII, are also divided in isotropy types as Kasner’s and the relation to the vanishing 
of the different sets of tensors carries over, in this case without exceptions. The (b.2), (b.l), (c), (a) types can be 
defined precisely setting to zero one further 4-tensor at a time; K = 0, 72 = K = 0, R = 72 = M = 0. 

In four dimensional Einstein gravity {N = 1 pure Lovelock) we can analyze much more general cosmological models. 
These have been classified long ago by Bianchi. Kasner metrics correspond to the Bianchi type I models in this 
classification, and represent the asymptotic solutions within those models. The situation is much more complicated 
for more general models, yet Kasner solutions still play a preeminent role. Bianchi type II models turn out to have 
Kasner asymptotics both to the past and the future. The Kasner solutions connected through a Bianchi type II 
trajectory define the so called Kasner map. For more general models, the BKL conjecture proposed that the Universe 
close to the initial singularity undergoes a series of oscillations, transitions between Kasner epochs where the expanding 
and contracting directions exchange their roles. A precise realization of this conjecture is given by the Mixmaster 
Universe. This represents the asymptotic behavior of Bianchi types VII and XI and is given by the iteration of the 
Kasner map, giving in this way the sequence of Kasner epochs. This sequence is infinite in four dimensions and the 
dynamics has been shown to be chaotic m- This chaotic behavior disappears in higher dimensional Einstein gravity. 

Non-trivial Kasner metrics in four dimensions correspond to type (b) solutions in our classification (type (b.l) are 
just flat Kasner in that case). It is reasonable to expect that a similar behavior to that of the Mixmaster attractor may 
also appear for other even dimensional pure Lovelock theories. It would be one more feature of d = 4 Einstein gravity 
respected by pure Lovelock theories in all dimensions. Our classification could also be regarded as the necessary first 
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step in such more general analysis. The obvious next step would be trying to generalize a result analogous to that 
of Bianchi type II in four dimensions and define a generalized Kasner map. For that we would have to introduce 
curvature in the spatial slices of our model (equivalently a nontrivial Lie group structure). The analogue of Bianchi 
II models would be a deformation in just a given 3-dimensional subspace; i.e. t^'dxi —> e® with i = 1,2 .. .d — 1 and 

[e^, e^] = n{t)e^ 

in the usual notation, the rest of the commutation relations being zero. This is the simplest possible modification of 
the Lie group structure. If a structure such as that of the Kasner map exists also in this case, this could also provide 
new examples of chaotic maps, interesting objects in their own right. The rich geometric structure we have found in 
the space of type (b) solutions would perhaps come in handy in performing such analysis. 


Appendix: Closed timelike curves (CTCs) from complex Kasner exponents 


The Kasner type metrics that are the subject of this paper are solutions of pure Lovelock gravity independently 
of the signature considered. The discussion of the next paragraphs will involve just the form of the metric and not 
the details of the equations, in that respect all that is said in this section applies also to Einstein gravity on any 
dimension. The choice of signature used in the previous sections corresponds to the description of the near region to 
a spacelike singularity. We can also describe timelike singularities in a similar fashion, we just have to exchange the 
roles of the time and one of the spacelike coordinates, t <->■ x. The form of the solutions, i.e. the possible values of 
the exponents is independent of the choice of signature, at least for real values of these exponents. 

Formally one can also choose complex values for the pi, as soon as they verify the relevant algebraic equations. The 
simplest such case is that of a complex conjugate pair of exponents, the rest being real. This is also an option in 4d 
Einstein gravity for which we can solve for two of the exponents in terms of the third as 


P± 


1 -P3 ± \/(l -P3)(3p3 + 1 ) 
2 


(A.l) 


These are real for a small range of values of ps € [—1/3,1], otherwise they have an imaginary part and we can write 
p± = a+bi, for real a, b. Notice that this is the range of values considered in the Lifshitz-Khalatnikov parametrization. 


Pi 


—u 

1 U + V? 


I -I- u 
1 u-\- v? 


u(I -I- u) 

1 + u + v? 


(A.2) 


with real parameter u. Even though u may take any real value the exponents always belong to the range pi S [—1/3,1]. 
Starting with Riemannian signature with a complex Kasner metric of the form. 


ds^ = e'^'^dr^ + e'^P'^dxj -f [e^^^-^dxl + e-^^^^dxl) , 

i^l,2 


(A.3) 


the term in brackets of the above metric is complex, this being the reason complex Kasner exponents have been 
disregarded in the past. We can however bring the metric to a real (Lorentzian) form performing a complex change 
of coordinates, xi = z = ti — it 2 ', X 2 = z = ti + it 2 . This yields. 


ds^ = e^^dr^ + ^ e^^''~dxj + 2 e^“’" [cos{2br){dt\ — dt^) + 2sin(25r) dti dt 2 ) (A.4) 

i/1,2 


even though we started off with Riemannian signature the presence of complex exponents already implies one timelike 
direction when the metric is written in real form. Notice however that the direction of the timelike direction in the 
(^ 1 ,^ 2 ) plane changes as we move along the radial direction r (see Figure]^ left). This makes a bit difficult to visualize 
the causal structure of this space and will determine some of its peculiar properties, namely the existence of closed 
timelike curves. Have we considered more than two complex exponents we would end up with more than one timelike 
direction. In some situations it is not even possible to bring the metric to any real form. 

For simplicity we will analyze the case a = I in detail. For o 7 ^ I some of the expressions will be a bit more 
complicated but the qualitative picture is still the same. We will consider just the {r,ti,t 2 ) part of the metric and 
neglect the conformal factor as it does not change the causal structure. Equivalently, we consider 

ds^ = dr^ + cos{2br){dt\ — dt\) + 2sin(2&r) dti dt 2 . 


(A.5) 
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We fix the future direction as that of positive t 2 at r = 0. This fixes the future direction for the rest of the space. It 
is orient able. For any r = const, plane it is easy to find the null lines as 


ti 


= 


' 1 =F sin(26r) 
1 ± sin(2&r) 


(A.6) 


These are ti = ±^2 for r = 0 and they rotate counterclockwise with constant opening as we increase r (see Figure]^. 
In addition to these, there is another family of null trajectories that is particularly simple to find. That is 


ti{r) = -cos{br) ; t 2 {r) =-sm{br) 


(A.7) 


that is future directed for f > 0. We can also have a future directed null trajectory going to lower values of r, i.e. 
r < 0, just changing the overall sign of ti, 

ti(r) = ^(2 - cos(6r)) ; t 2 {r) =-^ sm{br) (A.8) 

0 b 

where we added a constant such that both curves pass by (^, 0, 0). These two null curves coincide also for r = 2mT/b 
but they pass by these points in reversed order. The first curve is future directed from r = 0 to r = 27r/5 whereas the 
second is future directed from r = 27r/5 to r = 0. We can then take these two pieces and construct a closed future 
directed null curve (see Figure]^ right). This is the simplest case of such a curve, it is easy to work out many other 
examples. In the same way it is possible to demonstrate the existence of closed timelike curves, the null case is just 
some limiting case that is easier to analyze. In fact it is enough to add some small timelike pieces in between the 
trajectories going in both directions. This is true independently of the values of the complex exponents (equivalently 
a and b). 
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